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Qh! Abstract 

. We study the effects of Lorcntz non-invariance on the physical pion velocity at the critical 

temperature Tc in an effective theory of hidden local symmetry (HLS) with the "vector 
manifestation" fixed point. We match at a "matching scale" Am the axial- vector current 
' correlator in the HLS with the one in the operator product expansion for QCD, and present 

^ . the matching condition to determine the bare pion velocity. We find that the physical pion 

velocity, which is found to be one at T = Tc when starting from the Lorentz invariant bare 
HLS, remains close to one with the Lorentz non-invariance, U7r(Tc) = 0.83 — 0.99. This result 
is quite similar to the pion velocity in dense matter. 



1 Introduction 



Recent developments on effective field theory based on hidden local symmetry (HLS) suggest 
a state-of-the-art scenario for the chiral symmetry restoration at a large Nf, high temperature 
and/or high density, see Ref. [1] for a review. In the "vector manifestation (VM)" of the HLS 
theory, the p mesons become masslcss and the longitudinal components of the p mesons and 
the pions form a chiral multiplet at the restoration point [2]. It has been shown that the 
vector manifestation is realized for large Nf [2], at the critical temperature [3] as well as at the 
critical density [4] for chiral symmetry restoration. That the vector meson mass vanishes at the 
critical temperature/density supports the in-medium scaling of the vector meson proposed by 
Brown and Rho, "BR scaling" [5], and has qualitatively important influences on the properties 
of hadrons in medium. It has been shown [6] that the HLS with the VM predicts - and not 
posits - that the vector susceptibility xv equals the axial-vector susceptibility XA as required 
by chiral invariance and that the pion velocity Vj^ = 1 with the pion decay constants f^^O and 
^ as T — > Tg. This behavior differs drastically from the scenario predicted by the standard 
chiral theory [8] where the degrees of freedom at chiral restoration relevant to the iso-vector 
susceptibility are the pions after integrating out the scalar field in the 0(4) linear sigma model 
in the case of two-flavor QCD: The standard model prediction is that the pion velocity goes to 
zero at the chiral restoration point. 

One of the most important ingredients to realize the VM is the Wilsonian matching [7, 1] 
obtained from the following general ansatz: Integrating out quark and gluon degrees of freedom 
at a matching scale Am, we obtain the bare Lagrangian of the effective field theory (EFT). 
Then, we can determine the bare parameters of the bare Lagrangian by matching the EFT to 
the fundamental theory (QCD). Physical quantities are obtained by including quantum and 
thermal and/or dense effects through the renormalization group equations and thermal and/or 
dense loop. Since we integrate out the high energy modes, i.e., the quarks and gluons above 
A, in hot and/or dense matter, the bare parameters generically have the intrinsic temperature 
and/or density dependence [3, 4] which generally induces Lorentz symmetry violation in the bare 
EFT [6, 9]. 

In arriving at the HLS/VM results of Ref. [6], it was assumed that one can ignore Lorentz 
symmetry breaking in the bare HLS Lagrangian matched at a scale Am to QCD in medium. In 
this paper, we lift that assumption. This is made possible by a "non-renormalization theorem" 
recently proven by one of the present authors (C.S.) [10]. The observation is based on the 
existence of a new fixed point in the VM, namely that, the physical pion velocity does not receive 
any correction, either quantum or hadronic thermal, at the critical temperature. This means that 
it suffices to compute the pion velocity at the level of bare HLS Lagrangian at the matching 
point to arrive at the physical pion velocity at the chiral transition. It is important to study 
the physical pion velocity since this quantity is a dynamical object, which controls the pion 
propagation in medium through a dispersion relation. 

In the present work, we calculate the bare pion velocity by matching the axial-vector 
current correlators given by the HLS Lagrangian with Lorentz non-invariant terms taken into 
account, to the ones from the operator product expansion (OPE) in QCD. In Section 2, we briefly 
review HLS to define our notations and write down the HLS Lagrangian including the Lorentz 
symmetry violation. We present in Section 3 the Wilsonian matching condition to determine 
the bare pion velocity and its intrinsic temperature dependence in the low-temperature region. 
The bare pion velocity at the critical temperature is calculated in Section 4 by the Wilsonian 
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matching. Our results are summarized with a brief conclusion in Section 5. 



2 HLS Theory Without Lorentz Invariance 



The key observation of hidden local symmetry theory (HLS) is that any nonlinear sigma model 
defined in the coset space G/H is gauge-equivalent to a linear model possessing Ggiobai x -f^iocai 
symmetry. Here G = SU(A^j)j^ x SU(A'^^)j^ is the global chiral symmetry and H = SU(A'^^)y the 
HLS. In the HLS theory the basic ingredients arc the gauge bosons = p'^Ta of the HLS and 
two SU(A'^y)-matrix valued variables and ^r- They are parameterized as ^l,r = e'^'^^^" e^^'^^^^ , 
where tt = Tr^Ta denote the pseudoscalar Nambu-Goldstone (NG) bosons associated with the 
spontaneous breaking of G and a = a'^Ta the NG bosons absorbed into the HLS gauge bosons 
Pf^ which are identified with the vector mesons. and -Fq- are the relevant decay constants, and 
the parameter a is defined as a = F^/F!^. and transform as ^l,r(^c) ^ h{x)^i^^ji{x)gl^^, 
where h{x) G -ffiocai and gL,R & Ggiobai- The covariant derivatives of ^l,r are defined by D^^i, = 
dfj,^h — igPn^h + i(,L^^i, and similarly with replacement L ^ R, <-> 7^^, where g is the HLS 
gauge coupling, and and TZ^ denote the external gauge fields gauging the Ggiobai symmetry. 
The HLS Lagrangian is given by 



C = F^ti [a^^a>i]+F^tr 



+ ■^kin(/5/i) , 



where >Ckin(P(u) denotes the kinetic term of and 

«±,ii = (^/^^R-4T^^eL-ei)/(2i). 



(1) 



(2) 



As we have argued in the Introduction, the application of the Wilsonian matching in 
hot matter leads to the intrinsic temperature dependence of the bare parameters of the HLS 
Lagrangian which induces the Lorentz symmetry violation at the bare level [6, 9]. The Lorentz 
non-invariant HLS Lagrangian is constructed in Appendix A of Ref. [4]. At the lowest order, 
the Lagrangian is given by 

(^i^,bare)^«M^^ + ^^,bare-^^,bare idnu - U^U^) tr [a^d^] 
(-^(T,bare) '^/J.'^^i' + -^cr,bare-^iT,bare (ff^ii^ ~ U^jUy) tr 



+ 
+ 



avail 



5'L,bare 



UixUagvjB 



2g, 



r,bare 



tr 



(3) 



where F^ j^arg ,haxe) and -Fjj. 

bare) denote the bare parameters associated with the 
temporal and spatial decay constants of the pion (of the a). The unit four- vector needed to 
account for non-Lorentz-invariance in medium takes the value u'^ = (1,0) at the rest frame. 

Here it should be noticed that, due to the Lorentz symmetry violation, two variables and 
included in the 1-forms aj_ and ajf in Eq. (3) are parameterized as [10] 



Cl,r 



(4) 



*^ These scalars are not to be confused with the scalar that figures in two-fiavor hnear sigma model. 
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where ^^^.^ and ^^^^ are the bare parameters associated with the temporal decay constants 
of the pion and the a. 

The terms of 0{p'^) relevant to the present analysis are given by 



^Z2 — 



tr 



a/3 



(5) 



where the parameters ^2bare ^'^^ -^Jbare correspond in medium to the vacuum parameter 2;2,bare [7, 
1] at r = /X = 0. A'"' is defined by 



■^'"■=2 



1 r 



(6) 



where 7^''^ and L^"" are the field-strength tensors of the external gauge fields TZ^ and jC^: 



3 Matching Conditions for Bare Pion Velocity 

In this section, we present the matching conditions to determine the bare pion velocity including 
the effect of Lorentz symmetry breaking at the bare level. 

The Wilsonian matching is carried out by matching the vector and axial-vector current 
correlators derived from the HLS with those from the OPE in QCD. The axial-vector current 
correlator is defined by 

G7(go = i^n, q; T)6ab = fj^ dr j xe-'^^^-^+^-^^ {j^,{t, x)4,iO, 0))^ , (8) 

where J^^ is the axial- vector current, u;„ = 2mrT is the Matsubara frequency, (a, 6) = 1, ... , Nj — 
1 denotes the flavor index and ( the thermal average. The correlator for Minkowski momentum 
is obtained by the analytic continuation of qq. It is convenient to decompose the correlator into 

G^iqo, q) = q'PrG'iiqo, q) + q^PfOHqo, q) , (9) 
where we define q = \q\ and the polarization tensors are given by 




1' 



pf ^ -[r-^-^)-K'- (10) 

The bare Lagrangian is determined at the matching scale A through the matching and is ex- 
pressed in terms of the bare parameters as the sum of Eqs. (3) and (5). From this bare La- 
grangian, the current correlator at the matching scale is constructed as follows '^^i 

/ -\ 7r, bare TT, bare n L 

'^A(HLS)W0,9j — —To, 772 ^ ~ ^^2,bare ' 

\.% '^TT.barcy i 

*^The current correlator below A is obtained by including quantum and hadronic corrections through the 
renormalization group equations and thermal and/or dense loop. 
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GT I -\ TTjbare TTjbare r>^0'^2,bcire ^ ^2,bare (-\^\ 
A(HLS) Wo, gj ^2 ^ ^2 ' \^^) 

where T4-,bare = -^7? bare /-^iT bare t)are pion velocity. To perform the matching, we regard 

G^j^ as functions of —(^ and instead of go a^nd g, and expand G^j^ in a Taylor series around 
9 = 1^ = in /{.~Q^) as follows: 

Gi(-g^g-^) = Gr(-g^) + GT\-e)f + ■■■, 

Gli-q^f) = GT\-q') + G?^)(-g^)g-^ + • ■ ■ . (12) 

In the following, we determine the bare pion velocity V^,bare from and G^^^ via the 

matching. 

Expanding Q^^^^")^ Eq. (11) in terms of g^/(— g^), we obtain 

^(HLS)L(O)/ 2\ 7r,bare 7r,bare ^ L /i q\ 

) - ^^2,bare ' U-^j 

rpt rps /I T^2 \ 

G(HLS)L(1)/ 2\ 7r, bare TT, bare V-*- 7r,bare/ M A\ 

A {-Q ) {-q'^y ■ 

On the other hand, the correlator G'^ in the QCD sector to be given in OPE is more involved. 
Our strategy goes as follows. Since the effect of Lorentz non-invariance in medium has been 
more extensively studied in dense matter, we first examine the form of the relevant correlator in 
dense matter following Refs. [11, 12]. The current correlator G'^'^ constructed from the current 
defined by 

J^^^ = QlnQ^ or J^J = g757/xg, (15) 

is given by 



-coln\Q'\+j:^A-'- 



+ J2 ^[-g^^'q^^q^^ + g^'^'^q''q^'^ + q^'q^'^g''^'^ + gf^i^^ g''^'^Q'^] 

r=2fc = l 

22fc 

V«/*3 . . . „M2fc A2k+T,T 



T=2 fe=l ' ^ 



where = — g^. t = d — s denotes the twist, and s = 2A; is the number of spin indices of 
the operator of dimension d. Here A"'" represents the contribution from the Lorentz invariant 
operators such as A^'^ = ^ (^G^)^. ^^i^^-^'^j. and CfJ^^-Ji^,, arc the residual Wilson coefficient 
times matrix element of dimension d and twist r; e.g., ^)^f^2/i3M4 gi'^^n by 

<'i%M3M4 =i{ST{q^^,D^,Di_,,D^,q))^ , (17) 

where we have introduced the S3mibol ST which makes the operators symmetric and traceless 
with respect to the Lorentz indices. The general tensor structure of the matrix element of 
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^^i^-^'Jfc is given in Ref. [14]. For A; = 2, it takes the following form: 



A. 



a/3 Act 



PaPpPxPa - -jiPaPfigXa +PaPxgi3a + PaPa9xp + PpPxgaa 
+ P/3PagaX +PxPagap) + ^(5a/35ACT + 5aA5^CT + gaag0x)\M 



(18) 



For r = 2 with arbitrary k, we have [11]: 



A^ 

a 



2fc+2,2 
2k 

2fc+2,2 
2fe 



^2,2fc^2fc 'Ji-^'^J 



'^L,2fc^2fc + -^i-^' 



2,2fe^2fe 

rG aG 
L,2k^2k ) 



(19) 



where C^^k = ^ + 0{as), Cj] 



q,G 



2^ 0(05) and C^2fe ~ £^(05) (with the superscripts q and G 
standing respectively for quark and gluon) are the Wilson coefficients in the OPE [11]. The 
quantities A'^ and A^ are defined by 



2 / x" ^[q{x,fj,) + q{x,iJ.)] 



'■f 

Jo 



2 / dx x'^-^G{x,n) , 



(20) 



where q{x,fi) and G{x,iJ,) are quark and gluon distribution functions respectively. We observe 

that (16) consists of three classes of terms: One is independent of the background, i.e., density 
in this case, the second consists of scalar operators with various condensates {0)p and the third 
is made up of non-scalar operators whose matrix elements in dense matter could not be simply 
expressed in terms of various condensates {0)p. 

It is clear that Eq. (16) is a general expression that can be applied equally well to heat- 
bath systems. Thus we can simply transcribe (16) to the temperature case by replacing the 
condensates {0)p by {0)t and the quantities A^^^".];^^ and C^i^^'jfe by the corresponding quan- 
tities in heat bath. The higher the twist of operators becomes, the more these operators are 
suppressed since the dimensions of such operators become higher and the power of l/Q^ appear. 
Thus in the following, we restrict ourselves to contributions from the twist 2 (r = 2) operators. 
Then the temperature-dependent correlator can be written as 



Gr(go,g) = (gV-5'^V)-;^ 




a. 



1 + — In :V + 



9Q6 



6Q^ 



u,d,s 



+ [-g'^^q'^'q"^ + gl'l'^q'^q''^ + ql'qt'^ g^^l^^ + g'^l'^ g'^l'^Q'^] 



Lq4 ^1^2 



1 « 



(21) 



where G'^ is constructed from the axial-vector current associated with the iso-triplet channel 
defined by 



<^5m = 2^u'y5-^fj,u-d'y5-ffj,d), 



(22) 



5 



and we keep terms only up to the order of 1/Q^ for A^jt^/j^^,. The A" denote the SU{3) color 
matrices normalized as tr[A"A''] = 26°'^. Here we have dropped the terms with C^i^^2^^. in the 
non-scalar operators since they are of higher order in both 1 / {Q'^)"' and ag compared to the terms 
in the first line of Eq. (21). The temperature dependence of ^^f^^j ^^'^ ^^^^2/^3^4' implicit in 
Eq. (21), will be specified below. 

In order to effectuate the Wilsonian matching, we should in principle evaluate the con- 
densates and temperature-dependent matrix elements of the non-scalar operators in Eq.(21) at 
the given scale Am and temperature T in terms of QCD variables only. This can presumably be 
done on lattice. However no complete information is as yet available from model-independent 
QCD calculations. We are therefore compelled to resort to indirect methods and we adopt here 
an approach borrowed from QCD sum-rule calculations. 

Let us first evaluate the quantities that figure in Eq.(21) at low temperature. In low 
temperature regime, only the pions are expected to be thermally excited. In the dilute pion-gas 
approximation, {0)t is evaluated as 



{0)t ^ {O)o + 



3 

a=l-' 



2e(27r) 



-(7r»(p)|0|7r»(p))nB(e/T), 



(23) 



where e = ^Jp^^+m^ and is the Bose-Einstein distribution. As an example, we consider the 
operator of (r, s) = (2,4) that contributes to both G^^^ and G^^\ Noting that G\{qQ,q) = 
GAoo/^^ we evaluate Gaoo(?o,^)- 

(r=2,s=4)/ -\ _ 3 / d^P 16 Oa _0„/3 I _0 _a _0/3 I „0a _0^r)2i 



2e(27r)3 



where ^^^xa given by #^ 



(24) 



46,2(7r) _ 



P 



PaPfSPXPa - -^{PaPlSgXa + PaPxg^a + PaPaQxp + Pl3PX9o 



P 1 
+ PpPaQaX + PXPa9ap) + ^ {9af39Xa + gaX9pa + gaa9l3x)\ , (25) 

where carries the temperature dependence. Taking = 0, we see that the terms with 
and in Eq. (25) are zero. 

From Eqs. (21), (23) and (24), we obtain 

^(OPE)L(O). 2^ _ -i „z.^(OPE)(0) 




Q' 



*''For the general tensor structure of the matrix elements with a polarized spin-one target, say, along the beam 
direction in scattering process, see Ref. [16]. 
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4:7Tas 



u,d,s 



(26) 



^(OPE)L(i) ^g^j^gg following form 



G 



(OPE)L(l) 



32 „4^^ 4T(«+d) 

ids'" q^^^ ■ 



(27) 



We now proceed to estimate the pion velocity by matching to QCD. 
First we consider the matching between q(^^^)^(^) ^j^fj ^(OPE)l(o) j^qj^ ^-^^^ 
(26), we obtain 



d{- 



rpt rps 
^(HLS)L(O) _ -^7r,bare-^7r,bare 

) ^ - Q-2 ' 



d(-g2) 



G 



(OPE)L(O) 



-1 

8^ 



ri^^W^il^^ 3 1408 as{m)^T 
ttJ^ 3 g4 + ^ 27 



(28) 



Matching them at = A\j we obtain 



Tpt 

7r,bare 7r,bare 



1 

8^ 



a. X ^ 27r2 (^G^)r ^ ^3 1408 a. (^g)^ 



27 



15 Al, 



21 A^M ' 



= Go 



(29) 



Next we consider the matching between q^^^^^'I^W ^j^,^ G^2^^^^^^\ Prom Eqs. (14) and (27), 
we have 



TTjbare 7r,baj:e 



(1 ^,bare) 



32 4 Mu+d) 



(30) 



Noting that the right-hand-side of this expression is positive, we verify that 



^^,baj:e ^ 1 



(31) 



which is consistent with the causality. 

**The current correlator is expanded in positive power of in the HLS theory, while in negative power of 
in the OPE. Thus wc cannot match the whole Q'^-dependence between the HLS theory and OPE. However if 
there exists an overlapping area around a scale A, we can require the matching condition at that A. In fact, the 
Wilsonian matching at T = in throe flavor QCD was shown to give several predictions in remarkable agreement 
with experiments [7, 1]. This strongly suggests that there exists such an overlapping region. As discussed in 
Ref. [1], Am -C Ahls can be justified in the large Nc limit, where Ahls denotes the scale at which the HLS 
theory breaks down. Wo obtain the matching conditions in A^c = 3 by extrapolating the conditions in largo 
Nc- As we mentioned above, the success of the Wilsonian matching at T = with taking the matching scale as 
Am = 1.1 GeV shows that this extrapolation is valid. 
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The bare pion velocity can be obtained by dividing Eq. (30) with Eq. (29). What we 
obtain is the deviation from the speed of hght: 



S^^re = l-V^^^^^ = -—nj^A, . (32) 

This should be valid at low temperature. We note that the Lorentz non-invariance does not 

appear when wc consider the operator with s = 2, and that the operator with s = 4 generates 
the Lorentz non-invariance. This is consistent with the fact that including up to the operator 
with s = 2 is expressed as the function of only [15]. Equation (32) implies that the intrinsic 
temperature dependence starts from the 0{T^) contribution. On the other hand, the hadronic 
thermal correction to the pion velocity starts from the ©(T^) [9]. [There are O(T^) corrections 
to [/^]^ and [/^/^], but they are canceled with each other in the pion velocity.] Thus the hadronic 
thermal effect is dominant in low temperature region. At the critical temperature, the O(T^) 
corrections to the pion velocity are also cancelled as in the low-temperature region. Furthermore 
there are no 0{T^) corrections to cither [/^]^ or [/*/^]. Thus hadronic thermal corrections to 
the pion velocity are absent due to the protection by the VM. Therefore there remains only the 
intrinsic temperature dependence determined via the Wilsonian matching [10]. 

4 The Pion Velocity at Critical Temperature 

In this section, wc wish to evaluate the physical pion velocity at the critical temperature Tc 
starting from the Lorentz non- invariant bare Lagrangian. According to the non-renormalization 
theorem [10], the bare velocity so calculated should correspond to the physical pion velocity at 
the chiral transition point. Now, in the VM, bare parameters are determined by matching the 
HLS to QCD at the matching scale Am and at temperature T = Tc- 

We begin with a summary of the pion velocity found in the HLS/VM theory with Lorentz 
invar iance [6, 10, 9]. The pion velocity is given by [6]#^ 

"'^^^^ - F,2(0) + Reni(g-g-;r) ' ^''^ 

where IL^^{q, q;T) is the axial- vector two-point function [6] that represents hadronic thermal 
corrections. The two-point function Ilj^ is decomposed into four components, 

nf = u^u^Ili + {gf"" - u''u'')Ui + P^'^U^ + P^^Ul . (34) 

The VM dictates that if one ignores Lorentz non-invariance in the bare Lagrangian in medium, 
the pion velocity approaches the speed of light as T ^ Tc [6, 9]. 

In the following, we extend the matching condition valid at low temperature, Eq. (32), 
to near the critical temperature, and determine the bare pion velocity at Tc. As we discussed 
in the previous section, we should in principle evaluate the matrix elements in terms of QCD 
variables only in order for performing the Wilsonian matching, which is as yet unavailable from 
model-independent QCD calculations. Therefore, we make an estimation by extending the dilute 
gas approximation adopted in the QCD sum-rule analysis in the low-temperature region to the 
critical temperature with including all the light degrees of freedom expected in the VM. In 

*®Note that this definition is equivalent to the one used in Ref. [9] at one-loop order. 
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the HLS/VM theory, both the longitudinal and transverse p mesons become massless at the 
critical temperature since the HLS gauge coupling constant g vanishes. At the critical point, 
the longitudinal p meson which becomes the NG boson a couples to the vector current whereas 
the transverse p mesons decouple from the theory because of the vanishing g. Thus we assume 
that thermal fluctuations of the system are dominated near not only by the pions but also 
by the longitudinal p mesons. In evaluating the thermal matrix elements of the non-scalar 
operators in the OPE, we extend the thermal pion gas approximation employed in Ref. [14] 
to the longitudinal p mesons that figure in our approach. This is feasible since at the critical 
temperature, we expect the equality Al{Tc) = A1{Tc) to hold as the massless p meson is the 
chiral partner of the pion in the VM It should be noted that, although we use the dilute 
gas approximation, the treatment here is already beyond the low-temperature approximation 
adopted in Eq. (23) because the contribution from p meson is negligible in the low-temperature 
region. Since we treat the pion as a massless particle in the present analysis, it is reasonable to 
take Al{T) ~ AJ(r = 0). We therefore use 

A^iT) Al{T) ^ Al{T = 0) for T ~ T^. (35) 

The properties of the scalar operators giving rise to the condensates are fairly well un- 
derstood at chiral restoration. We know that the quark condensate must be zero at the critical 
temperature. Furthermore the value of the gluon condensate at the phase transition is known 
from lattice calculations to be roughly half of the one in the free space [13]. We therefore can 
use in what follows the following values at T = T^. 

{qq)T = 0, {—G^)t ~ 0.006GeV^ . (36) 

Including the contributions from both pions and massless p mesons, Eq. (27) can be expressed 
as 



^(OPE)L(l) ^ 32 ^4r' 



105 " Q' 

Therefore from Eq. (32), we obtain the deviation (^bare as 



(37) 



_ 1 32 , 

^TT.bare 



(38) 



This is the matching condition to be used for determining the value of the bare pion velocity 
near the critical temperature. 

To make cl rough. 6Stiniclte of (^bare? US6 

^^(n+d)^^ = 1 GeV) = 0.255 [14]. This value 
is arrived at by following Appendix B of [14] . ^dn^'^'* is defined by 

(7r|^7^,L>^,...L>^„g(/x)|7r) = (-z)"-'(Pmi -^'^n - traces) (/.) , (39) 

where 



A<i\p) = 2 /' dxx^-\q{x,p) + {-lTq{x,p)] . 
Jo 



(40) 



*®We observe from Refs. [19, 20] that (/^ = 2.4 GeV) ~ A^'-"~^'^\n = 2.4 GeV) even at zero temperature. 
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For any charge state of the pion (tt^, tt''", tt ), An^^'^^^ {n = 2,4) can be written in terms of the 
nth moment of valence quark distribution V^(/x) and sea quark distribution S'^{ix) [14], 

A^("+'^)(/x)=4C(/x) + 8S:(/x), (41) 

where 

= [\xx''-^v'^{x,ii), 
Jo 

SI = [ dxx'^-^s'^{x,ii). (42) 
io 

Simple parameterizations of the valence distribution v'"{x,n) and the sea distribution s^{x,ijl) 
can be found in Ref. [17]- see Ref. [18] for the updated results - where the parton distri- 
butions in the pions arc determined through the tt-N DreU-Yan and direct photon produc- 
tion processes. With the leading-order parton distribution functions given in [17], we obtain 
^^(u+d) _ g at = 1 GeV [14]. For the purpose of comparison with the lattice QCD 
result [19], we need to calculate the value at /x = 2.4 GeV; it comes out to be A^^^^^'^'^ = 0.18. 
The value ^4^"'^''^ = 0.18 is slightly bigger than the ~ 0.13 calculated by the lattice QCD [19], 
while it is a bit smaller than the ~ 0.22 of Ref. [20]. Note that the value ~ 0.18 agrees with 
the one determined by lattice QCD [19] within quoted errors. Using ^2^"+'^^ = 0.972 and 
^n{u+d) _ Q 255 a^Q(j fQj. ^]^g range of matching scale (Am = 0.8 — 1.1 GeV), that of QCD scale 
{Aqcd = 0.30 - 0.45 GeV) and critical temperature (Tc = 0.15 - 0.20 GeV), we get 

<Jbare(rc) = 0.0061 -0.29, (43) 

where the dependence of ^3 4'^'^^ is ignored as it is expected to be suppressed by more than 
1/A^. Thus we find the bare pion velocity to be close to the speed of light: 

K,barc(rc) =0.83-0.99. (44) 

Thanks to the non-renormalization theorem [10], i.e., v.,r(Tc) = V^,bare(rc)) we arrive at the 
physical pion velocity at chiral restoration: 

t,^(Tc) = 0.83-0.99. (45) 



5 Conclusion 

In this paper we lifted the assumption made without justification in the previous paper [6] 
about the ignorable role of Lorentz symmetry breaking in the bare Lagrangian which led us to 
the conclusion that the pion velocity at the chiral phase transition equals the speed of light. 
This feat of doing away with the assumption was made possible by the non-renormalization 
theorem [10] that states that in the HLS/VM theory, the VM protects the pion velocity from 
quantum as well as hadronic thermal corrections (at least at one-loop level) and hence the bare 
pion velocity obtained by matching to QCD at a matching scale Am remains un-renormalized 
by corrections at the critical temperature. By using information available from QCD sum 
rule calculations made in heat bath, we found that the pion velocity at the chiral transition 
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temperature remains close to the speed of light. This result is drastically different from the 
result obtained in the standard chiral theory [8] for the two-flavor QCD which predicts that 
the pion velocity should go to zero at the critical point. The crucial difference between the 
HLS/VM result and the standard chiral theory result lies in the degrees of freedom that figure 
at the phase transition: What accounts for the drastic difference in the prediction is the vector 
mesons becoming massless in the former at the VM fixed point to which the system is driven 
by temperature. 

It is interesting to compare the result obtained here in heat bath to a similar result 
obtained in dense medium, i.e., v^^ 1, as the system approaches the critical point. In [21], 
it was found that a dense matter is driven to chiral restoration by the change in the skyrmion 
background characterized by the expectation value in "sliding vacua" (SVEV in short) of the 
scalar field that represents the "soft glue" in QCD trace anomaly. It turns out that the SVEV of 
the scalar field relevant for this process is locked to the quark condensate (qq) in such a way that 
the melting of the quark condensate at chiral restoration precisely corresponds to the melting of 
the "soft glue" associated with part of the QCD trace anomaly. Remarkably the pion velocity is 
found to approach the speed of light in dense medium at the critical density in a way analogous 
to what is found in this paper for the same quantity at the critical temperature. This suggests 
that the vanishing SVEV of the scalar representing the soft glue in dense medium might be 
playing the role analogous to the VM fixed point in HLS/VM theory wherein the vector meson 
mass goes to zero. It remains to be seen whether the same result is obtained with HLS/VM in 
dense medium. 
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